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Abstract. To measure spin-dependent parton distribution functions in the production of bosons 
at the Relativistic Heavy Ion Collider, an accurate model for distributions of charged leptons from 
the W boson decay is needed. We present single-spin lepton-level cross sections of order ^{(X^) for 
this process, as well as resummed cross sections, which include multiple parton radiation effects. 
We also present a program RhicBos for the numerical analysis of single- and double-spin cross 
sections in y* , , and Z" boson production. 

The measurement of longitudinal spin asymmetries in the production of bosons 
at the Relativistic Heavy Ion Collider (RHIC) will provide an essential probe of spin- 
dependent quark distributions at high scales [1, 2]. At pp center-of-mass energy 
y/s = 500 GeV, about 1.3 x 10^ and W bosons will be produced by the time the 
integrated luminosity reaches 800 pb^^ . Due to the parity violation in the Wqq coupling, 
this process permits non- vanishing single-spin asymmetries A^{£,), defined here as 

dc>{p-^p-^WX) _ da{p^p-^WX) 

, da^p^'p^w -xj ^ Where ^=^,^,3;,,;?^;,.... (1) 

+ dl 

Here y^, is the rapidity of the W boson; and Pj^ are the rapidity and transverse 
momentum of the charged lepton from the W boson decay in the lab frame, respectively. 
The lowest-order expression for the asymmetry Aj^{y^,) with respect to the rapidity 
of the W boson is particularly simple if the absolute value of y^, is large. In that 
case, A^(jjy) reduces to the ratio Aq{x)/q{x) of the polarized and unpolarized parton 
distribution functions [3]. Furthermore, A^(^) tests the flavor dependence of quark 
polarizations. 

The original method for the measurement of W boson production at RHIC is based 
on the direct reconstruction of the asymmetry A^(jjy) [1]. Unfortunately, such recon- 
struction is obstructed by specifics of the detection ofW^ bosons at RHIC. First, RHIC 
detectors do not monitor energy balance in particle reactions. Due to the lack of informa- 
tion about the missing momentum carried by the neutrino, the determination of y^ is in 
general ambiguous and depends on assumptions about the dynamics of the process. Sec- 
ond, due to the correlation between the spins of the initial-state quarks and final-state 
leptons, the measured value of Aj^{y^,) is strongly sensitive to experimental cuts im- 
posed on the observed charged lepton. This feature is illustrated in Fig. 1, which shows 
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Figure 1. Dependence of the asymmetry Ajj-^) on the cuts imposed on the momentum of the observed 
antilepton in the boson production p^p {W^ lvi)X. The asymmetry is calculated using the 
resummation method described in the paper. The GRSV standard set [4] of the polarized PDF's was 
used. The error bars are calculated according to Eq. (13) in Ref. [1] assuming the integrated luminosity 



^ = 800 pb ' and beam polarization 70%. 



the asymmetry in the boson production calculated without constraints on yj and 
(solid line), and with constraints 1.2 < \yi\ < 2.4, pjj > 20 GeV (circles) and < 1, 
Pj,j > 20 GeV (boxes). According to the Figure, there is a substantial difference between 
the asymmetries calculated with and without selection cuts. This difference arises due 
to the different dependence of the unpolarized and polarized cross sections on angular 
distributions of the leptons in the W boson decay. For instance, at the lowest order in 
W+ production 

-Au{Xa)d{Xi^){l-COSe*)^ + Ad{Xa)u{Xi^){l+COSe*)^ 
dyw^yi)- u{Xa)d{Xf^) {I- COS 0*)^ + d{Xa)u{Xi^) {I + COS 0*)^ ' 

where x^^j^ = {Q/ y/s)e^yw^ and 0* is the polar angle of the antilepton in the rest frame 
of the W boson. Since is related to cos 0*, as 
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it is clear that restrictions on the range of integration of yj strongly affect A^{y^). 

Since the only straightforward signature of the W bosons at RHIC is the observation 
of secondary charged leptons, it is important to understand differential cross sections of 
spin-dependent W boson production at the lepton level. Given that the radiative correc- 
tions are sizeable (~ 30%) both in the numerator and denominator of Eq. (1), and that 



the measurement results will be used in the next-to-leading order (NLO) PDF analysis, 
it is necessary to derive these cross sections at NLO accuracy (i.e., at order ^(a^)). 

Furthermore, most of the W bosons are produced with small, but non-zero trans- 
verse momenta. Such non-zero qj. is acquired through radiation of soft and coUinear 
partons, which cannot be approximated by finite-order perturbative calculations. In or- 
der to obtain reliable predictions for differential cross sections, dominant logarithmic 
terms a^ln'" (^^ / Q^) (where <m<2n — I) associated with such radiation should be 
summed through all orders of the perturbative series. In our work [5], we performed a 
complete lepton-level study for the production of W"^, 7*, and ZP bosons for arbitrary 
longitudinal polarizations of incident protons. This study combined the ^(oCg) contri- 
butions with the all-order sum of small-^^ logarithmic corrections. The resummation of 
the logarithms ln"'(^^/2^) was performed with the help of the impact parameter space 
(^-space) resummation technique [6]. It extended the methodology developed for the 
unpolarized vector boson production [7] to the spin-dependent case. 

Our study goes beyond those in the previous publications [8, 9, 10] in several aspects. 
It presents the fully differential NLO cross section at the lepton level, which was not 
available before. The resummed single- and double-spin cross sections for the produc- 
tion of on-shell vector bosons were presented earlier in Ref . [10]. We have derived a more 
complete resummed cross section, which also accounts for the decays of vector bosons, 
i.e., for spin correlations in the final state. The lepton-level cross section includes several 
additional angular structure functions, which do not contribute at the level of on-shell 
vector bosons. Moreover, resummation is needed not only for the parity-conserving an- 
gular function 1 +cos^ 6*, which contributes to the on-shell cross section, but also for the 
parity- violating angular function 2 cos 0*, which affects angular distributions of the de- 
cay products. The estimate of parity-violating contributions is more complicated since 
it involves 7^ -matrices and Levi-Civita tensors both from the electroweak Lagrangian 
and spin-projection operators. As a result, special care is needed to treat finite terms 
that arise in the factorization of coUinear poles in d ^ 4 dimensions. We perform the 
calculation using the dimensional reduction and find that the resummed cross sections 
satisfy helicity conservation conditions for the incoming quarks. In addition, the coef- 
ficient functions in Ref. [10] were obtained in a non-conventional factorization scheme 
and cannot be used with the existing PDFs. In contrast, our results are fully consistent 
with the MS factorization scheme. 

The resummed cross sections are incorporated in a numerical program RhicBos for 
Monte-Carlo integration of the differential cross sections [11]. We are not able to discuss 
all aspects of our numerical study in this short report. However, as an example we show 
lepton-level asymmetries A^{yi) for various cuts on pj-i (Fig. 2). We find that these 
asymmetries can be accurately measured for both and bosons. These directly 
observed asymmetries can efficiently discriminate between different PDF sets; hence, 
they provide a viable alternative to the less accessible asymmetry Aj^{yyy). 

It is also useful to study distributions with respect to the transverse momentum 
of the charged lepton, not only because they are sensitive to the PDF's, but also be- 
cause they probe in detail dynamics of the QCD radiation. As was discussed above, 
the transverse momentum distributions for vector bosons are affected by the multiple 
parton radiation, which can be described only by means of all-order resummation. In 
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Figure 2. Asymmetries A^(y^) for various selection cuts on in (a) boson production and (b) 
boson production. The asymmetries are derived using the Gehrmann-Stiriing PDF sets A and B [12] and 
GRSV-2000 standard PDF set [4]. Statistical errors are estimated as in Fig. 1. 



addition, the distributions at very small qj- are sensitive to nonperturbative contributions 
characterized by large impact parameters b > I GeV^^. As a result, the shape of the 
lepton-level distribution do/dpj^ around its peak at about = M^jl is affected by 
both perturbative and nonperturbative QCD radiation. Remarkably, the shape of the Ja- 
cobian peak can be predicted by the theory, even though it cannot be calculated at any 
finite order of a^. The prediction is possible because the perturbative soft and coUinear 
contributions are systematically approximated in the resummation formalism. The non- 
perturbative contributions currently cannot be derived in a systematical way; but there 
is substantial indirect evidence (spin independence of the perturbative soft radiation, 
quark helicity conservation) that such contributions may be practically independent of 
the proton spin and type of the vector boson. The impact of the nonperturbative con- 
tributions is illustrated in Fig. 3, which shows the number of events for the difference 
da{p^p)/dpj.f —da{p^p)/dpj,j of single-spin cross sections at ^ = 800 pb^^. This 
rate was calculated using two parameterizations [13, 14] of the nonperturbative part, 
which were found in the unpolarized vector boson production. It can be seen that the sen- 
sitivity to the nonperturbative input is small, but, nonetheless, visible near the Jacobian 
peak. For comparison, we also included the ^{oCg) finite-order cross section calculated 
using the phase space slicing method. The finite-order curve substantially deviates from 
the resummed curves, and, moreover, its shape can be drastically modified by varying 
the phase space slicing parameter In contrast, the resummed curve is determined 
unambiguously once a parameterization of the nonperturbative input is obtained from 
the double-spin y* production, single-spin or double-spin Z'^ production, or even un- 
polarized W production. Needless to say, the resummation predictions for the shape of 
the Jacobian peak, which directly follow from fundamental principles of QCD, must be 
tested at RHIC. 
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Figure 3. The single-spin charged lepton transverse momentum distribution for boson production 
discussed in the main text. The nonperturbative parts of the resummed cross sections were calculated 
using the Ladinsky-Yuan parameterization [13] (solid) and the most recent Brock-Landry Gaussian 2 
parameterization [14] (dashed). The G'{a^) finite-order cross section is shown as a dotted line. 
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